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We examine the possibihty that a significant component of the energy density of the universe has an 
equation-of-state different from that of matter, radiation or cosmological constant (A). An example 
is a cosmic scalar field evolving in a potential, but our treatment is more general. Including this 
component alters cosmic evolution in a way that fits current observations well. Unlike A, it evolves 
dynamically and develops fluctuations, leaving a distinctive imprint on the microwave background 
anisotropy and mass power spectrum. 

PACS number(s); 98.80.-k,95.35.+d,98.70.Vc,98.65.Dx,98.80.Cq 

Inflationary cosmology predicts that the universe is 
spatially flat and that the total energy density of the uni- 
verse is equal to the critical density. This prediction is 
consistent with current measurements of the cosmic mi- 
crowave background (CMB) anisotropy and may be ver- 
ified with high precision in the next generation of CMB 
satellite experiments. At the same time, there is growing 
observational evidence that the total matter density of 
the universe is significantly less than the critical density.tl 
If this latter result holds and the CMB anisotropy estab- 
lishes that the universe is flat, then there must be an- 
other contribution to the energy density of the universe. 
One candidate that is often considered is a cosmological 
constant. A, or vacuum energy density. The vacuum den- 
sity is a spatially uniform, time- independent component. 
Cold dark matter models with a substantial cosmological 
constant (ACDM) are among the models which best fit 
existing observational data.EJ However, it should be em- 
phasized that the fit depends primarily on the fact that 
the models have low matter density and are spatially flat; 
the fit is not a sensitive test of whether the additional en- 
ergy contribution is vacuum energy. 

In this paper, we consider replacing A with a dynam- 
ical, time-dependent and spatially inhomogeneous com- 
ponent whose equation-of-state is different from baryons, 
neutrinos, dark matter, or radiation. The equation-of- 
state of the new component, denoted as w, is the ratio 
of its pressure to its energy density. This fifth contri- 
bution to the cosmic energy density, referred to here as 
"quintessence" or Q-component, is broadly defined, al- 
lowing a spectrum of possibilities including an equation- 
of-state which is constant, uniformly evolving or oscil- 
latory. Examples of a Q-component are fundamental 
fields (scalar, vector, or tensor) or macroscopic objects, 
such as a network of light, tangled cosmic strings .0 The 
analysis in the present paper applies to any component 
whose hydrodynamic properties can be mimicked by a 
scalar field evolving in a potential which couples to mat- 
ter only through gravitation. In particular, we focus on 
equations-of-state with — 1 < w < because this range 
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fits current cosmological observations 
has motivated several investigational' t 
with w < in which a spatially uniform distribution has 
been assumed, e.g., a decaying A or smooth component. 

In this Letter, we begin by arguing that a smoothly 
distributed, time-varying component is unphysical — it 
violates the equivalence principle. Hence, predictions of 
CMB and mass power spectra which have not included 
fiuctuations in the new component are not valid. We 
outline the general conditions needed to have w < con- 
sistent with the equivalence principle and stable against 
catastrophic gravitational collapse. We show that an 
evolving scalar field automatically satisfies these condi- 
tions. We then compute the CMB and mass power spec- 
tra for a wide, representative class of models (see Figure 
1 and 2a) , taking careful note of the effects of the fluctua- 
tions in the Q-component. We show that the fluctuations 
leave a distinctive signature that enables a Q-component 
to be distinguished from dark matter and cosmological 
constant and makes it possible to resolve its equation- 
of-state. Comparing to observations of the CMB and 
large-scale structure (Figure 2b) , we identify a new spec- 
trum of plausible models that will be targets for future 
experiments. 

Introducing a dynamical energy component is at least 
as well motivated by fundamental physics as introducing 
a cosmological constant. In fact, the theoretical preju- 
dice based on fundamental physics is that A is precisely 
zero; if it is non-zero, there is no conceivable mechanism 
to explain why the vacuum density should be compara- 
ble to the present matter density, other than arguments 
based on the anthropic principle. On the other hand, 
dynamical fields abound in quantum gravity, supergrav- 
ity and superstring models {e.g., hidden sector fields, 
moduli, pseudo-Nambu-Goldstone bosons), and it may 
even be possible to utilize the interaction of these fields 
with matter to find a natural explanation why the Q- 
component and matter have comparable energy densities 

As noted above, a number of studiesi'&i have as- 
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sumed a "smooth" (spatiahy uniform), time-dependent 
component with arbitrary equation-of-state (sometimes 
called xCDM) which does not respond to the inhomo- 
geneities in the dark matter and baryon-photon-neutrino 
fluid. In computing the CMB anistropy, one finds a near- 
degeneracyEl with ACDM for a wide range of w. In this 
Letter, we show that the degeneracy is significantly bro- 
ken when fluctuations in the Q-component are included. 
However, it is important to realize that the smooth and 
Q-scenarios are not competing models. A smooth, time- 
evolving component is ill-defined, since the smoothness 
is gauge-dependent, and unphysical, because it violates 
the equivalence principle to ignore the response of the 
new component to the inhomogeneities in the surround- 
ing cosmological fluid. Hence, a fluctuating, inhomoge- 
ncous component is the only valid way of introducing an 
additional energy component. 

There have been some discussions in the literature of 
an energy component consisting of a dynamical, fl." 
ating cosmic scalar field evolving in a potentialjj'li 
or an energy conTOaricjjit evolving according to a specific 
equation-of-state.Qla't3 In the case of a scalar field, the 
CMB anisotropy and power spectcum were computed by 
Coble et al for a cosine potentialM-and by Ferreira and 
Joyce for an exponential potential.Hil Here we go beyond 
these isolated examples to explore the range of possibili- 
ties and the range of imprints on the CMB and large-scale 
structure. 

The class of cosmological models we consider in this 
work are spatially flat, Friedmann- Robertson- Walker 
(FRW) space-times which contain baryons, neutrinos, ra- 
diation, cold dark matter and the Q-component (QCDM 
models). The space-time metric is given by ds^ = 
a?{rj){—drf' -\- dx^) where a is the expansion scale fac- 
tor and T] is the conformal time. For the purposes of 
CMB and mass power spectrum prediction, we model the 
Q-component of the fluid as a scalar fleld, Q, with self- 
interactions determined by a potential V{Q). In an ideal 
adiabatic fluid with w < 0, a concern would be that the 
sound speed is imaginary, = w < 0, and small wava^ 
length perturbations are hydrodynamically unstable.t3 
However, a scalar field is not an ideal fluid in this sense. 
The sound speed is a function of wavelength, increasing 
from negative values for long wavelength perturbations 
and approaching =^-.1 at small wavelengths (smaller 
than the horizon size) Jlil Consequently, for the cases con- 
sidered here, small wavelength modes remain stable. In 
general, our treatment relies only on properties of w and 
Cj, and so it applies both to cosmic scalar fields and to 
any other forms of matter-energy (e.^., non-scalar fields 
or topological defects) with w and which can also be 
obtained by a scalar field and some potential. Otherwise, 
we make no special assumptions about the microscopic 
composition of the Q-component. 

We implicitly assume that any couplings to other fields 
are negligibly small, so that the scalar field interacts with 



other matter only gravitationally. The average energy 
density is pq = ^Q'^ + V and the average pressure is 
Pq ~ 2a^Q'^ ^ y 1 where the prime represents d/drj. As 
Q evolves down its potential, the ratio of kinetic {^^Q'^) 
to potential {V) energy can change; this would lead to 
a time- varying equation-of-state, w = pq/ pq. In this 
paper, we consider constant and time- varying w models 
for which w = pq/ pq € [—1,0], since this range includes 
models which best fit current observations. 

We have modified numerical Boltzmann codes for com- 
puting CMB anisotropy and the mass power spectrum. 
In one approach, we use an explicitly- defined equation- 
of-state {'w{rf)), which spans all models but requires that 
the equation-of-state as a function of time be known. All 
properties of the model are determined by the proscribed 
wirf) = w\a{ri)\ and 



PQ = ^qPo exp 
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where Po and are respectively the critical energy den- 
sity and scale factor today. An effective scalar potential 
which produces this pq is computed as a byproduct. A 
second approach uses an implicitly-defined equation-of- 
state in which we specify the potential V{Q) directly, 
which is useful for exploring specific models motivated 
by particle physics. 

A series of codes have been developed for synchronous 
and for conformaLNewlonian gauge by modifying stan- 
dard algorithmsJl3'E3EZl_In the synchronous gauge, the 
line element is given byliS 



ds^ ~ a^{ri)[—dri^ + {'fij -\- hij)dx''dx^] 
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where 7^ is the unperturbed spatial metric, and hij is 
the metric perturbation. The scalar field fluctuation SQ 
obeys the equation 



SQ" + 2-6Q' - V^6Q + a^Vqql 
a 



(3) 



Here, h is the trace of the spatial metric perturbation, as 
described by Ma and Bertschinger.E3 

We have found the observable fluctuation spectrum to 
be insensitive to a broad range of initial conditions, in- 
cluding the case in which the amplitudes of SQ, 6Q' were 
set by inflation. All of the examples shown in this paper 
have SQ = SQ' — initially as measured in synchronous 
gauge. With this choice of initial conditions, the imprint 
of 6Q on the CMB anisotropy would be negligible if we 
did not include the response of Q to the matter density 
perturbations through the h'Q' source term in Eq. (H). If 
Q' — 0, as occurs for ACDM, the scalar field remains un- 
perturbed. In the implicit approach where V is given, the 
energy density, pressure, and momentum perturbations 
are 



;Q'6Q' + VqSQ, Spq = -^Q'6Q' - Vq6Q, 
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These quantities must be included in the evolution of 
h, the total fluctuation density contrast S, and the total 
velocity perturbation v. For cases where w{r]) is provided 
explicitly, the same relations may be applied mhere an 
effective V is computed by the code from w(?7)Jla 

Although we have examined many models ,113 we re- 
strict ourselves here to a representative spectrum. Figuxe 
1 illustrates CMB studies showing COBE-normalizedO 
CMB anisotropy spectra assuming adiabatic initial con- 
ditions in the matter with tilt n = 1, Hubble con- 
stant Ho = 100ft km s-i Mpc"i with h = 0.65, and 
r2f,/i^ = 0.02, where fib is the ratio of the baryon density 
to the critical density. Because the models are spatially 
flat, the angular size of the sound horizon at last scat- 
tering is nearly the same for all our models and, conse- 
quently, the acoustic oscillation peaks are_aLnearly the 
same fs as in standard CDM (SCDM).B0 However, 
the shape of the plateau for low £ and the shapes of the 
acoustic peaks at high £ are distorted by three effects: 
a corrJaination of early and late integrated Sachs- Wolfe 
effectsEil due to having a fluid component with w ^ 
and the direct contribution of Spq on the CMB. Fig- 
ures l(a)-(b) show how the spectra vary with fig, the 
ratio of the energy density in Q to the critical density, for 
flxed w. The behavior is different for different equations- 
of-state. For w = —1/6 the flrst acoustic peak rises and 
then steadily decreases as fig increases; for w ~ —2/3, 
the acoustic peak rises uniformly as fig increases. Fig- 
ure 1(c) compares predictions assuming a smooth, non- 
fluctuating component (a case we argued is unphysical) 
with the full spectrum including fluctuations in Q, illus- 
trating the substantial difference at large angular scales. 
Figure 1(d) shows spectra for fixed fig — 0.6 and varying 
w. When fluctuations in Q are not included, thje curves 
for — 1 < w < — 1/2 are very nearly degenerate,El but the 
degeneracy is substantially broken in Figure 1(d) by the 
contribution of Q fluctuations at large angular scales at 
a level much greater than cosmic variance. Figure 1(e) 
illustrates rpaults for a time- varying w{ri) obtained using 
exponentialEj and cosined scalar potentials. The expo- 
nential potentials are of the form V{Q) = m'^ exp[—PQ] 
where Q — Q' — initially; for the examples shown, 
(m, (3) = (0.00358 eV, 11.66 m-i) for the case with 
w{tjo) = -1/6 and (m, /3) = (0.00243 eV, 8.0 m"!) for 
10(770) = —2/3 as expressed in Planck mass (rrip) units. 
For the cosine potential, V{Q) = m'^{l+cos[Q/ f]), where 
(to, /) = (0.004B5eV, 0.1544mp) with Q = 1.6/ and 
Q' = initiallyE Two of the curves are especially in- 
teresting because they compare the CMB anisotropy for 
two models with the same value of w today but with 
different values of w in the past. One model has con- 
stant w = —1/6 whereas the other model, based on a 
scalar field rolling down an exponential potential, has w 



evolving from -1 to -1/6. The two models produce eas- 
(4)ily distinguished CMB power spectra even though the 
equation-of-state is the same today due to the difference 
in past evolutionary history and the direct contributions 
of SQ. Hence, the evolution of w can also be determined 
from CMB measurements. Also shown are examples of 
other exponential and cosine potentials compared with 
SCDM which illustrate that the CMB anisotropy is sen- 
sitive to V{Q). Figure 1(f) illustrates an example of CMB 
polarization, which has similar variations with parame- 
ters as the temperature anisotropy. 

Figure 2 illustrates the mass power spectrum predic- 
tions. The key feature to note in Figure 2(a) is that 
the power spectrum is sensitive to all parameters: the 
value of w; the time-dependence of w; the effective po- 
tential V{Q) and initial conditions; and the value of Qq. 
Hence, combined with the CMB anisotropy, the power 
spectrum provides a powerful test of QCDM and its pa- 
rameters. An important effect is the suppression of the 
mass power spectrum and erg (the rms mass fluctuation at 
8h^^ Mpc) compared to SCDM, which makes for a bet- 
ter fit to current observations .0 The grey swath in Figure 
2(b) represents the, constraint from x-ray cluster abun- 
dance for ACDM;E3 in general, the constraint is weakly 
model-dependent. 

While QCDM and ACDM both compare well to cur- 
rent observations of CMB and of large-scale structure 
today, QCDM has advantages in fitting constraints from 
high red shift supernovae, gravitational lensing, and 
structure formation at large red shift (z « 5). Con- 
straints based on classical cosmological tests on A (wa — 
— 1), such as supernovae and lensing, are significantly re- 
laxed for QCDM with w « —1/2 or greater .Ej Another 
property of QCDM (or ACDM) is that structure growth 
and evolution ceases when the Q-component (or A) be- 
gins to dominate over the matter density. Comparing 
QCDM and ACDM models with Qq = JIa, this cessa- 
tion of growth occurs earlier in QCDM. For larger values 
of w, the growth ceases earlier. Hence, more large scale 
structure and quasar formation at large red shift are pre- 
dicted, in better accord with deep red shift images. 

In conclusion, we find that the "quintessence" hypoth- 
esis fits all current observations and results in an imprint 
on the CMB anisotropy and mass power spectrum that 
should be detectable in near-future experiments. Its dis- 
covery could indicate the existence of new, fundamental 
fields with profound implications for particle physics, as 
well as cosmology. A number of follow-up studies are 
underway including: a quantitative analysis to deter- 
mine how well future CMB experiments can resolve a 
Q-component and how well one can simultaneously re- 
solve other cosmic parameters; numerical simulations of 
galaxy formation and evolution in the presence of a Q- 
component; and search for quintessential candidates in 
models of fundamental physics. 

We thank P.J.E. Peebles, J.P. Ostriker and M. White 



3 



for useful comments. This research was supported by the 
Department of Energy at the University of Pennsylvania, 
DE-FG02-95ER40893. 



[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 
[9] 
[10] 

[11] 

[12] 

[13] 

[14] 
[15] 

[16] 
[17] 
[18] 
[19] 
[20] 



[21] 

[22] 
[23] 



See, for example, J. P. Ostriker and P.J. Steinhardt, Na- 
ture 377, 600 (1995), and references therein. 
David Spergel and Ue-Li Pen, Astrophys. J. 491, 67 
(1997). 

P. J. Steinhardt, in Critical Dialogues in Cosmology, ed. 
N. Turok, (World Scientific, 1997). 

Michael S. Turner and Martin White, Phys. Rev. _D56, 
R4439 (1997). 

K. Coble, S. Dodelson, and J. A. Frieman, Phys. Rev. 
i?55, 1851 (1997). 

J. Frieman, C. Hill, A. Stebbins, and I. Waga, Phys. Rev. 
Lett. 75, 2077 (1995) 

L. Wang, R. R. Caldwell, J.P. Ostriker, and P.J. Stein- 
hardt, in preparation. 

V. Silveira and I. Waga, Phys. Rev. D56, 4625 (1997). 
C. Wetterich, Astron. Astroph. 301, 321 (1995). 

B. Ratra and P. J. E. Peebles, Phys. Rev. D37, 3406 
(1988). 

Pedro G. Ferreira and Michael Joyce, Phys. Rev. Lett. 
79, 4740 (1997). 

Takeshi Chiba, Naoshi Sugiyama, and Takashi Naka- 
mura, Mon. Not. R. Astron. Soc. 289, 5 (1997). 
See, for example, J. C. Fabris and J. Martin, Phys. Rev. 
D55, 5205 (1997). 

L. Grishchuk, Phys. Rev. D50, 7154 (1994). 

R. Crittenden, J.R. Bond, R.L. Davis., G. Efstathiou, 

and P.J. Steinhardt, Phys. Rev. Lett.71, 324 (1993). 

C. -P. Ma and E. Bertschinger, ^p. J. 455, 7 (19 95). 



U. Seljak and M. Zaldarriaga, a,stro-ph/9603033, 
R.R. Caldwell, R. Dave, P.J. Steinhardt, in preparation. 
E.F. Bunn and M. White, Ap. J.480, 6 (1997). 
For reviews of CMB anisotropy power spectrum physics 
and properties, see P. J. Steinhardt, in Proceedings of 
the Snowmass Workshop on the Future of Particle Astro- 
physics and Cosmology, ed. by E.W. Kolb and R. Peccei, 
(World Scientific, 1995); IJMPA AlO, 1091-1124 (1995). 
W. Hu, in The Universe at High-z, Large Scale Struc- 
ture and the Cosmic Microwave Background, eds. E. 
Martinez-Gonzalez and J.L Sanz (Springer Verlag); 
Wayne Hu, Naoshi Sugiyama, and Joseph Silk, Nature 
386, 37 (1997). 

P.T.P. Viana and A.R. Fiddle, Mon. Not. Roy. Astron. 

Soc. 281, 323 (1996). 

P.J. Steinhardt, Nature 382, 768 (1996). 



o 



8 


(a)w = -1/6 


(b)w='-2/3 'A 




SCDM 


n„=o.8 /a 


6 






4 










£2„=o.2 M . 


2 




SCDM 




1 


1 1 1 


8 


(c)w = -1/3, n„=o.6 


(d) % = 0.6 
ACDM 


6 


— smooth 

component /A 


W = -2/3 $ W 
W = -1/2 ///^ 


4 


Q-component / \ 


W = -1/3 I'/ \| /• A ■ 


2 




SCDM^^ yl^fc- 






1 — 1 — 


8 


1 — 1 

(e) time-varying w(ti) 


(f) = 0.6, Polarization 


6 


e"PQ ;w„ = -M . 
— e-PQ;w. = -,/6.^;,^^ 


(XIOO) )i 
ACDM O'l) 


4 


-.-cos(Q/f);w.= . . 
w = -1/6 V/v^f- 




2 


SCDM f 11 


n 







10 



100 



1000 



100 1000 10 

Multipole moment (/) 

FIG. 1. CMB power spectrum, £{£ -\- l)Ci/2TT vs. £ where 
Ce is the multipole moment, illustrating that the spectrum 
changes significantly as a function of w and Qq. Thick 
solid lines represent SCDM; thick dashed lines correspond to 
ACDM. For thin solid lines, legends list the parameters in se- 
quence according to the height of the curve beginning from the 
topmost. The variation with fig is shown in (a) and (b); (c) 
compares the predictions for a Q-component with fluctuations 
{SQ 7^ 0) to a smooth component (5Q=0). The differences 
due to 5Q at large angular scales change COBE-normalization 
which is responsible for the substantial variation of the acous- 
tic peaks with w, as shown in (d). Panel (e) shows the results 
for time- varying w{ri) as determined by specifying a potential 
V{Q) and initial conditions, where V{Q) is an exponential or 
cosine functional of Q in the examples shown. Included also is 
the prediction for a constant w — —1/6, which differs from the 
results where w{ri) is time-varying and io(?7o) ~ —1/6. In (f) 
is shown the polarization (amplified 100- fold) for w = —1/2 
compared to SCDM and ACDM. 
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FIG. 2. (a) Variation of mass power spectrum for some 
representative QCDM examples, (b) The variation of erg with 
rig. For ACDM, rig is ^Ia- The suppression of erg in QCDM 
compared to SCDM makes for a better fit with current ob- 
servations. The grey swath illustrates constraints from x-ray 
cluster abundance. 
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